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The author and colleagues have produced a collection of 60 test problems which offer a 
realistic perfornriance test of the currently available automatic codes for eigenvalues of the 
classical Sturm-Liouville problem. We describe a Fortran implementation and the 
considerations that went into its design. A novel feature is that (almost) all the code 
defining one problem is textually contiguous in the Fortran text, unlike for example the 
DETEST package for ODE initial-value solvers where the def ... 
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High order domain decomposition methods using a basis of Legendre polynomials, known 
variously as ''spectral elements" or ''p-type finite elements/' have become very popular. 
Recent studies suggest that accuracy and efficiency can be improved by replacing 
Legendre polynomials by prolate spheroidal wave functions of zeroth order. In this article, 
we explain the practicalities of computing all the numbers needed to switch bases: the 
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Divided differences are enormously useful in developing stable and accurate numerical 
formulas. For example, programs to compute ^(x)-f(/) as might occur in integration, can 
be notoriously inaccurate. Such problems can be cured by approaching these 
computations through divided difference formulations. This paper provides a guide to 
divided difference theory and practice, with a special eye toward the needs of computer 
algebra systems that should be programmed to deal ... 
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In the last 20 years there has been nearly a six-order-of-magnitude increase in computing 
power, and it is predicted that another five-to-slx-orders-of-magnitude increase will occur 
by the year 2000 [1]. Although supercomputers such as the Cray-1 and Cyber 205 can 
perform certain vector operations with burst rates of up to 100 million floating-point 
operations per second (Mflops) [2], it is remarkable that a popular and affordable 
personal computer is only three-orders-of-magnitude slower. ... 
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Publisher: ACM Press 

Full text available- ■P ^pdff63146 KB) Additional Information: full citation , abstract, references, citings, ind^ 
^ terms 

Many geometric problems like generalized Voronoi diagrams, medial axis computations 
and boundary evaluation involve computation and manipulation of non-linear algebraic 
primitives like curves and surfaces. The algorithms designed for these problems make 
decisions based on signs of geometric predicates or on the roots of polynomials 
characterizing the problem. The reliability of the algorithm depends on the accurate 
evaluation of these signs and roots. In this paper, we present a {\em naive ... 
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Publisher: ACM Press 

Full text available: pdf(17.07 MB) Additional Information: fuil citation , abstract 

Level set methods, an important class of partial differential equation (PDE) methods, 
define dynamic surfaces implicitly as the level set (iso-surface) of a sampled, evolving nD 
function. The course begins with preparatory material that Introduces the concept of using 
partial differential equations to solve problems in computer graphics, geometric modeling 
and computer vision. This will Include the structure and behavior of several different types 
of differential equations, e.g. the level set eq ... 
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A probability density on a finite-dimensional Euclidean space Is orthounimodal with a 
given mode if within each orthant (quadrant) defined by the mode, the density is a 
monotone function of each of its arguments individually. Up to a linear transformation, 
most of the commonly used random vectors possess orthounimodal densities. To generate 
a random vector from a given orthounimodal density, several general-purpose algorithms 
are presented; and an experimental performance evaluation illustr ... 
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